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Abstract 

In previous publications we have analyzed the strong and electromagnetic decays of 
heavy mesons and heavy baryons in a formalism which incorporates heavy-quark and chiral 
symmetries. There are two possible symmetry-breaking effects on the chiral dynamics of 
heavy hadrons: the finite- mass effects from light quarks and the 1 / rriQ corrections from heavy 
quarks. In the present paper, chiral-symmetry-breaking effects are studied and applications 
to various strong and radiative decays of heavy hadrons are illustrated. SU(3) violations 
induced by chiral loops in the radiative decays of charmed mesons and charmed baryons 
are compared with those predicted by the constituent quark model. In particular, available 
data for D* decays favor values of the parameters in chiral perturbation theory which give 
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predictions for D* decays close to the quark model results except for the D^'^. Imphcations 
are discussed. 
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1. Introduction 



The dominant decay modes of many heavy hadrons, which contain a heavy quark, are 
strong decays with one soft pion emission and/or electromagnetic decays. This is a con- 
sequence of the heavy quark symmetry [1,2] of QCD: mass differences among the different 
spin multiplets of the ground-state heavy mesons and heavy baryons are generally small. 
An ideal framework for studying the low-energy dynamics of heavy hadrons is provided by 
the formalism in which the heavy quark symmetry and the chiral symmetry are synthesized 
[3-10]. However, symmetry considerations alone in general do not give any quantitative pre- 
dictions unless further assumptions are made. Fortunately, all the unknown parameters in 
the Lagrangian depend only on the light quarks and are calculable from the nonrelativistic 
quark model. In Refs.[3,9,10] (for later convenience, Ref.[3] and Ref.[10] will be denoted as 
paper I and paper II respectively, henceforth), we have explored in detail the predictions of 
this theoretical formalism on strong decays, radiative decays, and heavy-flavor-conserving 
nonleptonic decays. 

In this and a preceding paper, we would like to examine various symmetry-breaking 
corrections to the strong and radiative decays of heavy mesons and baryons. There are two 
different kinds of symmetry-breaking effects on the chiral dynamics of heavy hadrons: the 
finite-mass effects from the light quarks and the l/mg corrections from the heavy quarks. 
In paper II, we have already incorporated one of the l/mg effects, namely the magnetic 
moment of the heavy quark, into the formalism for describing the electromagnetic (Ml) 
decay of heavy hadrons. This is because the charmed quark is not particularly heavy, and 
hence the contribution due to its magnetic moment cannot be safely neglected. 

There are two strong motivations, among others, for promoting a systematic study of 
both the l/mg corrections and the effects of chiral symmetry breaking. First, we have 
calculated in paper II the decay rates of D* D^. When combined with our prediction for 

the strong decays D* Dtt given in paper I, we arc able to predict the branching ratios 
for the D* decays. Agreement is excellent between theory and the most recent experiment 
of CLEO II [11]. Nevertheless, our predicted total width for D*+ is rtot(^*+) = 150 
keV,* which is to be compared with the upper limit rtot(-D*"'") < 131 keV published by the 

7*r The difference between this number (see Table I) and the result 141 keV obtained in Ref.[10] is due to 
the fact that here we have used a more accurate pion momentum. 
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ACCMOR Collaboration [12]. We thus urge the experimentalists to perform more precision 
measurments of the D* total width. Therefore, it becomes urgent to analyze the symmetry- 
breaking effects on the strong decays D* — > Dn and the radiative decays D* — > D^. This is 
particularly so should the aforementioned upper limit for Ftotl^*^) be confirmed by future 
experiments. Secondly, to the lowest order chiral Lagrangian, there exist several SU(3) 
relations among the radiative decay amplitudes of heavy mesons and heavy baryons, for 
example, A{D*^ D^l) = ^{D'^'^ D^l)- Therefore, observation of different rates for 
the decays D*^ D'^'j and D*'^ D~^l (after taking into account mass differences between 
and Dt^*^) will clearly signal the SU(3) flavor symmetry breaking due to the hght 
current quark masses. In paper II, the magnetic moments of heavy hadrons are related to 
the magnetic moments of the constituent quarks in the nonrelativistic quark model. In this 
approach, SU(3) violations in radiative (Ml) decays arise from the constituent quark mass 
differences. In the present work, we shall follow the orthodox approach of chiral perturbation 
theory to treat SU(3)-breaking effects: The coupling constants in the Lagrangian are treated 
to be invariant and all SU(3) violations of interest are induced by chiral loops. A comparison 
of these two different approaches of SU(3) breaking is given in the end. Existing data for the 
D* decays, though limited, favor a set of parameters such that chiral perturbation theory 
and the quark model give similar results for the D* decays, but not for D*+ Dfj. In 
addition, even if the upper limit for the total width of D*^ [12] is to be taken seriously, it 
does not pose a difficulty in the chiral perturbation theory approach. For the heavy baryons, 
lack of any data prevents us from making a meaningful search for an optimum choice of 
parameters in chiral perturbation theory. 

Since both symmetry breaking effects require a careful and thorough study, we focus on 
the SU(3) symmetry breaking in this paper. A detailed investigation of l/mg corrections 
has been carried out in a preceding paper [13]. Schematically, the general effective chiral 
Lagrangian in chiral perturbation theory (ChPT) involving heavy hadrons has the chiral 
expansion^ 

Ci + C2 + Cs + ---, (1.1) 

t For pure Goldstone-boson fields, the chiral expansion of chiral Lagrangians has the familiar form 

£ = £2 + A + -Ce + • • • . 
In this case, higher order chiral terms are suppressed by powers of P^/A^ or m^/A^. 
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where the subscript denotes the sum of the number of derivatives acting on the Goldstone 
fields and the velocity-dependent heavy hadron fields, and the power of Goldstone-boson 
mass squared in the Lagrangian. (Recall that the lowest-order chiral-symmetry-breaking 
terms are linear in light quark masses and hence quadratic in Goldstone-boson masses.) The 
higher-order chiral Lagrangians are suppressed by powers of p/A^ or m^/A^, where p (rri) 
is the momentum (mass) of the Goldstone bosons, and A,^ is a chiral symmetry breaking 
scale ~ 1 GeV [14]. Therefore, perturbation theory makes sense if p and m are not too large 
compared to A^^. Chiral corrections of interest usually receive two contributions: one from 
the chiral loops generated by Ci, and the other from the higher order tree Lagrangian term 
£2- Loop contributions can be either finite or divergent. The divergent part typically has 
the form (see Appendix A) 

--7,+ln47r + l + ^ln^ (1.2) 

in the dimensional regularization scheme, where e = 4 — n, and A is an arbitrary renormal- 
ization scale. In ChPT, all divergences from chiral loops induced by Ci^ will be absorbed 

into the counterterms which have the same structure as that of C2 [15]. Denoting the bare 
parameters of C2 by /i, the renormalized parameters fl are then given by [16] 

/[(A) = + ^ - 7. + ln47r + 1) A^ (1.3) 

with 7j being calculable coefficients. We see that although the lowest-order chiral Lagrangian 
C\ is scale independent, renormalized higher-order effective Lagrangians do depend on the 
choice of A, reflecting the non-renormahzability nature of ChPT. Of course, physical ampli- 
tudes should be independent of the renormalization scale; that is, the A dependence from 
chiral- loop is exactly compensated by the A dependence of local counterterms in £2- As 
the renormalized parameters fl are unknown and must be determined from experiment, 
we will thus concentrate in the present paper the chiral corrections due to meson loops. 
Furthermore, we will choose A ~ A^ to get numerical estimates of chiral-loop effects. 

Chiral loop corrections to some heavy meson processes have been discussed by other 
authors [8, 17-20]. In the heavy baryon sector, the only chiral loop effects studied so far are 

X The quadratic divergence, if exists, usually amounts to renormalizing Ci. In the dimensional regular- 
ization scheme we only focus on those logarithmic divergences. 
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those in Ref.[17] on the semileptonic decays. In the present paper, we attempt to make a 
systematic and complete investigation of chiral-loop corrections to the strong and electro- 
magnetic decays of heavy hadrons. They arise from the chiral-symmetry-breaking terms in 
(1.1) and hence vanish in chiral limit. The leading chiral-loop effects we found are nonana- 



light quark mass). Furthermore, they amount to finite- light-quark-mass corrections to the 
coupling constant, say g, in Ci. Schematically, 



where ^^(A) is the relevant renormalized coupling constant in C2; that is, the chiral-loop 
and £2 contributions have the same structure as the g term in Ci except that they vanish in 
chiral limit. This point will be elaborated on again in Section 2. In this work we shall only 
keep the nonanalytic loop effects. 

The present paper is organized as follows. In Section 2 we calculate chiral corrections 
to the strong decay P* — > Pn and the radiative decay P* — > P7, where P* and P refer to 
1~ and 0~ ground-state vector and pseudoscalar mesons, respectively. As noted earlier, the 
l/friQ effect due to the magnetic moment of the heavy quark is included. Similar chiral-loop 
calculations are presented for heavy baryons in Section 3, except there we adopt velocity- 
dependent "superfields" which combine spin-| and spin-| sextet baryon fields together [6,17]. 
Computation becomes much simplified in this compact notation. A by-product of our in- 
vestigation of Sections 2 and 3 is the confirmation at the one loop level of an exact QCD 
result (see [13] for example) that the coupling constants due to heavy quarks in the Ml 
transitions of both heavy mesons and baryons are not modified by the light quark dynamics. 
In Section 4 we consider applications of our results to the strong and electromagnetic decays 
of charmed mesons and charmed baryons. SU(3) violation induced by chiral loops for the 
radiative decays is compared with that predicted by the nonrelativistic quark model. Section 
5 contains conclusions and outlook. Two Appendices are devoted to some technical details. 



lytic in the forms of m/ K-^ and (m^/A^) ln(A^/m^) (or m'q and ruq Inm^, with ruq being the 




(1.4) 



6 



2. SU(3)-Symmetry-BrecLking Corrections to the Chiral Dynamics of Heavy 
Mesons 



In this section we shall study the chiral-symmetry-breaking effects on the strong and 
electromagnetic decays of heavy mesons. To begin with, we recall the lowest-order gauge 
invariant chiral Lagrangians (1.2.20) and (II. 2. 19) for heavy mesons 

C^^^pp^ = -2iMp,P{v)v ■ DP\v) + 2iMp.P*''{v)v ■ DP*\v) 

+ AM^P{v)P\v) + f^jMpMp* [P{v)A>'P*^\v) + P*{v)A''P\v)] (2.1) 
+ 2igMp. e^,xnP*^{vVA^P*''\v), 

and 

6^1p^ = yjMpMp.e^,^pv'^P'^P{v)[\d{i^Qi + ^Q^^) + d'Q!]F^-P\v) + h.c. 

1 (2.2) 

+ id Mp,F^,P''^{v)Yi^ --{E}QE,+E,Qi^)]P'^^\v), 



with AM2 = Ml, - Ml, 



p*t — n p*t _ n p*t 
p* _ n p* _ Tj p* 



where P and P* denote the ground-state 0^ and 1~ heavy mesons respectively which contain 
a heavy quark Q and a light antiquark g, and A^^ are the respective chiral vector and 
axial fields (see paper II for more detail), Q = diag(|, — is the charge matrix for the 

hght d and s quarks, Q' (or en) is the charge of the heavy quark, ^ = exp(-^M-) with 
the unrenormalized decay constant /o to be determined later, and M is the meson matrix 
of Goldstone boson fields 



M=> 



\/2 

TT' 



V K- K -JinJ 



(2.4) 



where A's are the Gell-Mann matrices normahzed by tr(A'*A^) = 25°'^. It should be stressed 
that the Lagrangians (2.1) and (2.2) are expressed in terms of velocity-dependent heavy 
meson fields. 
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The covariant derivative in Eqs.(2.1-2.3) is defined by 

D^P = d^P + V;P + ieA^{Q;P -PQ), 

(2.5) 

D^Pl ^df^Pl + V^pt _ ieA^{plQ! - QPl)^ 

with = ^[.^"I^D^^ + ^{D^j^^)'^]. The covariant derivative in contains the photon field A^ 
to incorporate the electromagnetic interactions of the Goldstone bosons; explicitly 

d^,i + ieA^[Q, i]. (2.6) 

In Eq.(2.2), v is the velocity of heavy mesons, the two couphng constants dQ and d"Q can 
be related to the magnetic moments of light constituent quarks in the quark model. The 
universal coupling constants d and d" are independent of the heavy quark masses and species. 
We have also included the d'Q! and 7Q' terms to account for the corrections due to the heavy 
quark masses when rriQ ^ 00. The coupling constants / and g in (2.1), d and d" in (2.2) are 
related by heavy quark symmetry as [3,10] 

/ = 2g, d" = -2d. (2.7) 

Moreover, the couplings d' and d"^ are fixed by heavy quark symmetry to be [10] 

d'^-w^, d"i^—. (2.8) 
2mQ ruQ 

We shall see that from the lowest-order Lagrangians given by (2.1) and (2.2), there is an 
SU(3) prediction: A{D*+ D+^) = A(P)*+ D+^). Therefore, observation of different 
rates for the radiative decays of D*'^ and D*+ (after taking into account the mass difference 
of and Dt^*^) will evidently signal the 5'[/(3)-breaking effects induced by the hght 

quark masses. In chiral perturbation theory the masses of light quarks are treated as a small 
perturbation and SU(3) violation in radiative decays is induced by chiral loops. In paper 
II, the unknown coupling constants d and d" in (2.2) are derived from the non-relativistic 
quark model and they are related to the masses of the constituent u, d and s quarks. As a 
consequence, SU(3) violation is already incorporated into the "effective" couplings d and d" 
in paper II. In the present work, we shall adopt the orthodox approach of chiral perturbation 

theory for treating SU(3)-breaking effects. As will be shown, leading chiral corrections have 

1 

nonanalytic dependence on niq of the form m| or niq Inm^. In Sec. 4 the above two different 
approaches for SU(3) violation in radiative decays of heavy hadrons will be compared. 
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The chiral symmetry breaking terms are given by 

>CcsB=f Sotr(A^tE + EtM) 

+aiMpP{v){^M^^ + ^^M^^)pHv) + a2Mp.P*{v){^M^C + C^M^^)P*^{v) (2-9) 
+a3MpP{v)P^{v)tr{M^E + E^M) + a4Mp*P*''{v)P*^ {v)tr{M^E + Il^M), 

where S = i?o = — 4 (qq) //q characterizes the spontantous breaking of chiral symmetry, 
and is a hght quark mass matrix 



M 



fmu \ 

ma 
\ msJ 



(2.10) 



The unknown dimensionless parameters ai arc expected to be of order unity (see Sec. 4 for 
more discussion). Each symmetry breaking term in (2.9) transforms either as (3,3) or (3,3) 
under SU{3)jj x SU{3)p. As the symmetry breaking effect is purely induced by the light 
quark masses, £cSB should be independent of the heavy quark mass and spin. By comparing 
(2.9) with (2.1), this immediately gives ai — —a.2 and 03 = — a4. As one can see from (2.9), 
>CcsB does not break the heavy-quark-symmetry (HQS) relations given by Eq.(2.7). In fact, 
chiral symmetry breaking should preserve HQS as light meson interactions have nothing to 
do with the heavy quarks. 

We shall first calculate the chiral-loop contributions to the strong decay P* Pj + tt" as 
depicted in Fig. 1, where Pi denotes {Qu, Qd, Qs) and likewise for P^ . The tree amplitude 
can be read out from (2.1) to be 

A[Pt ^ P,-7r«(g)]tree = MpMp.{\%j{e* ■ q). (2.11) 

Jo 

As discussed in the Introduction, the loop-induced logarithmic divergence is absorbed into 
the chiral-symmetry-breaking counterterms 

aP{v){iM^i + i'^Mi'^)A^P;\v) + bP{v)A^{iM^i + i^Mi'^)P;\v) + h.c. (2.12) 

Obviously, in the above expression, we should take ^ = 1 in {^M^^ -\- ^"^M^^) in order to 
describe the decay P* — > Ptt. As a consequence, the higher order contribution due to (2.12) 
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has the same structure as that of (2.11) except that the former vanishes in chiral hmit. 
As stressed in the Introduction, we will not consider higher-order Lagrangian effects in the 
present paper because of the unknown parameters a and b in (2.12). The sum of all chiral 
loop contributions shown in Fig. 1 gives rise to the effective coupling constant 

VZ2{P)Z2{P*)Z2{7:-) MPh/Z^ 

9es = g = g , (2.13) 

where Zi and Z2 are the vertex and wave-function renormalization constants, respectively, 
and the HQS relation Z2{P*) — Z2{P) has been applied. For simplicity, wave-function 
renormalization and mass counterterms of the heavy meson and Goldstone boson are not 
explicitly shown in Fig. 1, but the usual mass renormalization procedure is to be understood. 
To evaluate the Z2 renormalization constant, we note that the self-energy amplitude of, for 
example Pi, has the expression 

-lU{k) = -WJ2^P l^i^ ^ )ii / ToZ^ 772 2 , • NT /, , I • 1 ' ^^-^^f 

2/0 ^ J {^'^r{l^-m^^^+ie)[v-{l + k)+ie\ 

where k is the residual momentum defined by P^ — rnQV^ -\- k^, and m^b is the mass of the 
Goldstone boson TT^ in the loop. (For notation consistency, we will follow papers I and II to 
denote the momenta of the Goldstone boson and photon by q and k, respectively.) With the 
help of Eq. (A4) and the relation 

U{k) ^6m'^ -2{Z^^ -l){vk)Mp, (2.15) 

we obtain 

Sm^Pi) = Y.^X'^X%^Mp„ (2.16) 
-■-"^ f, JO 

and 

Z2(fl)=l + ^E(^'A.^l„^, (2.17) 

where A is an arbitrary renormalization scale. As discussed in the Introduction, divergences 
from chiral loops are absorbed into the unrenormalized parameters of higher-derivative chiral 
Lagrangians which are not written down here. It is straightforward to check that the HQS 
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relation Z2{P*) — Z2{Pi) holds. Note that the SU(3) invariant masses Mp and Mp* in 
Eq.(2.11) are unrenormalized masses; they are connected to the physical masses through the 
relations: 

M^ySPi) -Ml- 2aiMpMn - 2asMptrM + Sm\Pi), 

(2.18) 

M^^JP*) = Ml, - 2aiMp,Mii - 2a3Mp,tTM + Sm^{P*). 

It follows from Eqs. (2.9), (2.16) and (2.18) that the mass splitting of, say P3 = (Qs) and 
Pi — {Qu), is given by 

Likewise, for heavy vector mesons 

phvs plivs phvs phys 

Mp, -Mp, = Mp^ -Mp^ . (2.196) 

Consequently, a meaurement of the heavy-meson mass differences will provide information 
on the parameter ai. We will come back to this in Sec. 4. 

From the strong-interaction chiral Lagrangian for the Goldstone bosons 

C = iltr + ^tr (m^J: + ATst) (2.20) 

with A4 — diag(m^, m^, 2m^ — m^), the wave-function renormalization constants for the 
light pseudoscalar mesons are found to be 

yrz^^l-^(2e^ + e,), 

V^M =1 - lil^. + + le,), (2.21) 

^^=1-^(36,), 



with 



1 a2 

1 m2 , , 

"^In^, (2.22) 



327r2 /, 



"^i- 

1 ^^^A^_ 



" 327r2 /2 -2 



m 



V 



As to the Feynman diagrams Figs, la-le, the amplitude of the vertex diagrams Figs, la 
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and lb reads 



A{la + lb) = JL_^MpMp. J2(Xbx-X%-^ln^{s* ■ g), (2.23) 

h •'U Tr" 



where we have apphed Eq.(A3). To evaluate the seagull graphs of Figs. Ic-ld, we note that 
the Feynman rules for the vertices PPir'^Tr^ and P*P*tt°'tt^ are obtained from the kinematic 
terms in Eq.(2.1) together with Eqs. (2.3-2.6) and 

= -^(Md^M - d^MM) + • • • , (2.24) 
4/0 

where M is the meson matrix given by Eq.(2.4). The relevant Feynman rules are 



' Mp(A«A^ - A^A> • (g« - %) (2.25) 



for the vertex PPn^n^, and 

' Mp{X''X^-X^X^)v{qa-qb){e*-e'*) (2.26) 



4/o^ 

for P*P*7r°'7r^. Therefore, both vertices are proportional to v ■ {qa — qh)- As a consequence, 
seagull diagrams Figs. Ic and Id vanish since the required e* ■ q expression cannot be 
generated from the seagull amplitudes in which the linear q dependence is always in the v ■ q 
form. Beyond the heavy quark limit, seagull graphs do contribute but they are suppressed by 
factors of 1/Mp* relative to the tree amplitude. For example, to order 1/Mp*, the following 
p*p*^a^b vertex 



' (A«A^ - A^A«)Mp[(5« - %) • e'*{v ■ e*) - {qa - qb) ■ e*{v ■ e'*)] 



8/^ 



are no longer vanishing as v • £* 7^ and v ■ e'* ^ 

The 5-point vertex in Fig. le is generated from the Lagrangian (2.1) by expanding 
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to the third power of M 

A^, = d^M + — i=-3 [{d^M)M^ - 2M{d^M)M + M%M] +■■■. (2.27) 

A straightforward calculation yields 

r i(2e^ + e^), iorP*^Pj + n, 

A{le) = j-y/MpMp.{X%j{e* ■ q) x I i(f + + |e^), for ^ P,- + K, (2.28) 

[i(3e,), for/^*^P, + r^. 

Since the vertex corrections are lumped into the renormalization constant Zi by 

A{la + lb+lc+ld+ le) = (Zf 1 - 1)A{P* ^ Pj7r")tree, (2.29) 

it is clear from Eqs.(2.13), (2.21) and (2.28) that the contribution due to the 5-point ver- 
tex diagram Fig. le is exactly compensated by the wave-function renormalization of the 
Goldstone bosons. As a consequence, we have effectively 



with 



b 

Using the relations 



^i(la + l&)=l-fe^^^^e.- (2.31) 



J](A''A^)n/(7r'') = J2i>^'X%2f{n') = 3/(7r) + 2f{K) + \f{r,), 



b b 

b 



(2.32) 



it follows from Eq.(2.17) that the individual wave function renormalization constants for 

heavy mesons are given by 



Z2(Pi) = Z2{P2) =l + g^ Qe, + 3e, + ^e^^ , 
^2(^3) = 1 + 2^2(36^ + e^). 



(2.33) 



In Sec. 4 the above results will be applied to the strong decay D* ^ D + n. There we shall 
see that Z2 plays an essential role and the pion contribution is not negligible in spite of its 
small mass. 
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Thus far we have expressed all the results in terms of the unrenormalized decay constant 
/o- It can be related to the physical decay constants through the relations [16]: 



2ml+ml . ml 



A = /o ( 1 + 2e, + + 4 " , ^ L\ + 4-f L5 

/o H 

1 + + ^e, + 4 5!i±I=iLS + 2 f 1 . (2.34) 



2m^+m^ 4m^2ms + m 
1^ ''^'^11^ 



where rh = (m,, + mfi)/2, and the counterterm contributions denoted by the renormalized 
coupling constants L\ and (see Ref.[16] for notation) are included in Eq.(2.34). Note 
that the one-loop logarithmic corrections to the decay constants, the pion wave-function 
renormalization and the 5-point vertex diagram all share the same structure as they come 
from the same chiral loop diagram. The couplings L\ and L5 of the four-derivative chiral 
Lagrangian are dependent of the renormalization scale A. However, physical decay constants 
can be verified to be independent of the choice of A, as they should be [16]. Using the 
empirical values 

L^(A = m^) = (0 ± 0.5) X 10-^ L^(A = m^) = (2.2 ± 0.5) x lO'^ (2.35) 
obtained in Ref.[16] and the experimental value — 93 MeV [21], we find 

/o ~ 86 MeV. (2.36) 

We make a digression here to comment on the 1 / itlq effects. In the heavy quark effective 
theory, ^jraq corrections can be systematically studied by including the following higher- 
dimensional operators 

^""^ (2.37) 
2mQ 2 

where h^^^ is a velocity-dependent heavy quark field. For the case of P* — > Ptt decays, 
there are two different kinds of l/mq corrections. One arises from the consideration of the 
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subleading decay amplitude induced by the operators Oi and O2 



As,.[P*iv,e*)^P{v')7r%q)]^ j-q^{P{v')\zT J + O2(x)M^(0) |P*(^;, £*)) . 

(2.38) 

As shown in the preceding paper [13], the couphng constant g receives a correction of order 
Kcd/^q; that is, 

'A. 



5'eff = 9 



1 + 



'-QCD 



(2.39) 



This l/mg effect in chiral hmit is elaborated on in more details in Ref.[13]. The other 
kind of 1 / rriQ corrections comes from the combination of chiral-symmetry and heavy-fiavor- 
symmetry breakings, e.g., the seagull diagrams. Schematically, it is expected to be order 
of 

1 mK m^b 

7T^-k^^ 2.40 
{Any rriQ 

where the factor l/(47r)^ is associated with the loop- momentum integration, and m^b is the 
mass of the loop meson tt^. As this SU(3)-violating l/mg effect is very small, we will not 
pursue it further. 

We next turn to the electromagnetic decay of heavy mesons. One-loop contributions 
to P* ^ P + ^ are shown in Fig. 2. The Feynman diagrams with a photon coupled to 
the external P* or P line via charge coupling are vanishing, as expected since the reaction 
under consideration involves a magnetic Ml transition. Therefore, they are not displayed 
in Fig. 2. Although chiral-loop effects in the radiative decays have been recently discussed 
by Amundson et al. [8], only Fig. 2e was analyzed by them. Before proceeding to compute 
loop corrections, we note that the tree amplitude obtained from (2.2) is 

A[P*{v, £*) ^Pi + 7(^, s)] = Pie^uapki'e^v^e*^, (2.41) 

with 

Pi = -2y/ MpMp*{eid + cqd!) = e^MpMp*{eil3 + eq^'), (2.42) 

and Si — Qii- Though the unknown parameter (5 = —2d/e is not fixed by heavy quark 
symmetry, the other parameter /3' = —2d'/e is determined in the heavy quark effective 
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theory by the dimension-five operator 



with the result [7,8,10] 



O3 = ^hi'^\~eeQa^,F^^'')hi'^\ (2.43) 



P' = — . (2.44) 



It should be stressed that it is important to include this l/mg effect due to the magnetic 
moment of the heavy quark. As we have seen in paper II, the charmed quark contribution 
is significant and largely destructive in the radiative decays of D*^ and D^'^. In the non- 
relativistic quark model P is related to the constituent quark mass, namely P — 1/Mq [10]. 
Since in the spirit of chiral perturbation theory SU(3) violation is induced by chiral-symmetry 
breaking terms, it is interesting to compare the quark model results with those of the chiral 
Lagrangian approach. To the lowest order, it is evident that p2 = pz or 

A(P| ^ P27) = A{P^ ^ P37)- (2.45) 

As will be seen, this SU(3) relation is badly broken by chiral-loop effects. 
Referring to Fig. 2 the effective couplings become 

(ft)dl = P^f^ + *ft(2« + 26 + 2e), 

(2.46) 

where the subscript i refers to the process P* — > Pj7 and Zi accounts for the vertex renor- 
mahzation due to Figs. 2c-2d and 2f. We begin with Fig. 2e. In order to evaluate its 
amphtude, all the charged meson loops are exphcitly shown in Fig. 3. By virtue of Eq.(A2), 
we find* 



SPi{2e) 



ei V Stt 871 
1 f g"^ rriT, 



62 \8n /o 
1 f g'^ rrij. 



,2 



63 VStT 



■k In Ref.[8], /o is replaced by the physical decay constant for the pion loop and for the kaon loop. 
We will not make such a replacement in our calculation. 
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which are in agreement with Ref.[8]. As for the seagull diagrams Figs. 2a and 2b, we note 
that the electromagnetic four-point interactions of heavy mesons are described by 



JO JO (2.48) 

X {P{V)[Q, M]P;\v) + P;{v)[Q, M]P\v)), 



where use has been made of 



1 26 



(2.49) 



in (2.1) and (2.2). It is clear that the P*Pt^°'^ vertex is independent of the photon's momen- 
tum k. As a consequence, Fig. 2b with a photon emitted from the vertex on the r.h.s. does 
not contribute since no linear /c-dependent terms can be constructed from the numerator or 
the denominator after the loop mementa of the tt^ and P* are assigned to be / and p + I, 
respectively. Since the vertex P*P*7r"7 is proportional to e'^^'^'^v^SpS^e'*^ in the heavy quark 
limit, it is easily seen that the other seagull diagram Fig. 2a also vanishes. Its subleading 
contribution is of order (gf^//g)(l/MpN.) and hence negligible. 

We next focus on vertex corrections. With the magnetic P*P^ and P*P*^ vertices 
determined from (2.2), Figs. 2c and 2d lead to 

Zi(i^*^P,7;2c+2d)= 1-^ j;^^^§^6,., (2.50) 



and 



Z[{P: ^ Po; 2c + 2d)=l + ^(A^A^),,e^, = Z2(Pi), (2.51) 



where we have applied (2.17) and the relations (2.7) and (2.8). The contribution of Fig. 2c 
due to charge P*P*7 couphng is of order [g^ /?,2T:'^)[rn?^i / fQ){m^b / Mp) and hence neghgible. 
Note that because Z[{P* — > Pi^y) is exactly compensated by Z2{Pi), the parameter (3' of 
(2.46) does not get renormalized; this is a realization in chiral perturbation theory of the 
exact QCD result discussed in Ref.[13]. 
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It remains to compute the amplitude of Fig. 2f. The relevant electromagnetic 5-point 
vertex obtained from (2.2) is 



2/0^ 



^MpMp. e^.apk^'e^v^e^f^iM, [Q, M]]. 



(2.52) 



It is easily seen that only charged meson loops contribute; the result is 



Z,{p: Pa; 2/) = 1 + i $][A^ [Q, A^]],,g- 



(2.53) 



Working out (2.50) and (2.53) exphcitly, we find the total vertex corrections to be: 



Pit, 2c + 2d + 2/) = 1 + — 

ei 

P27; 2c + 2d + 2/) =1 + - 
P37; 2c + 2d + 2f) =1 + - 



(2.54) 



Putting everything together, we obtain the following effective couphngs for P* — > P^: 



pi = ey/MpMp* 
P2 = ey/MpMp* 



1/5 



Z2{Pl) 



eg 



3" Zi(P* ^ P17) mg Stt /2 



8^ U 



Z2{P2 



+ 



+ 



3"Zi(P|^P27) Stt /2 



P3 ^e^/MpMp* 



Z2{P3) 



^ eg ^ g'^ 



3"Zi(P|^P37) Stt 



(2.55) 



In the SU(3) limit, the light quark contributions to p are in the ratio 2 : — 1 : — 1 for the 
electromagnetic decays of Pj*, P|, and P|. Evidently, this relation is violated by the wave 
function and vertex renormalizations. 



3. SU(3)-Symmetry-BrecLk;ing Corrections to the Chiral Dynamics of Heavy 
Bciryons 

In this section we discuss the corrections to the strong and radiative decays of heavy 
baryons containing a heavy quark Q and two light quarks. The two light quarks form either 
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a symmetric sextet 6 or an antisymmetric antitriplet 3 in flavor SU(3) space. We will denote 
these spin ^ baryons as Bq and respectively, and the spin | baryon by Bq. Explicitly, 
the baryon matrices read as in paper I 



B, 



( 

1 



Q 







" 2 



\ 



/ 



(3.1) 



( 

1 

72^Q 



V2 Q 



\ J_'^/2 1 t:/" 



(3.2) 



^/2-Q ' 

and a matrix for Bq similar to Bq. The superscript in (3.1) and (3.2) refers to the value of 
the isospin quantum number 13. 

To perform chiral-loop calculations, we find that for sextet heavy baryons it is very 
convenient to work with the "superfields" which combine spin-| and spin-| sextet baryon 
fields into a single field [6,17] 



^^=^r-^(7'^+^'^)75i?6, 

S''-Bl^ + ^BQ^f>{l^ + v^l 



(3.3) 



where Bq^ is the Rarita-Schwinger vector-spinor field as appropriate for spin | particles. 
Feynman rules in terms of superfields become much simpler and there are fewer Feynman 
diagrams to evaluate. The most general gauge-invariant chiral Lagrangian for heavy baryons 
given by Eq.(II.3.8) reads 



'(1) 



- itT{S^'v ■ DS^) + -iiiTv ■ DT) + Amtr(5^5^) 
+ i\ gle^,afiMS''v''^^^S'') - V3^2tr(>SM^T) + h.c. 



(3.4) 



in superfield notation, where T = B^, Am is the mass splitting between the sextet and 
antitriplet baryon multiplets (for simplicity, we will drop the unknown Am term in ensuing 
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loop calculations), and 



D^B = df,B + V^S + SVj + ieQ'A^B + ieA^{Q, B} (3.5) 

for B — Si,, T. We wish to stress that the Lagrangian (3.4) is expressed in terms of velocity 
dependent baryon fields. Because of heavy quark symmetry, there are only two independent 
coupling constants in the low-energy interactions between the Goldstone bosons and the 
heavy baryons. Moreover, the nonrelativistic quark model predicts that [3]: 



91 = :^9: 92^ -y ^9 (3.6) 

with g being the axial vector coupling constant of a constituent quark, which is also the 
coupling constant appearing in the meson Lagrangian (2.1). Feynman rules can be easily 
derived from (3.4). Especially, the S and T propagators are simply given by i{—9^u + 
Vf^Vi,) /{v ■ k — Am) and i/v ■ k, respectively, where k is the residual momentum of the heavy 
hadron. The lowest order chiral-symmetry-breaking terms now have the form 



CcsB = Ai tr[^^(eA^^e + ^^M^^)S^] + A2 tiiS'^S^MM^i: + ^^M) 



(3.7) 



We begin by first examining the chiral-loop effects on the strong decays Bq — > Bq + tt'^. 
The tree amphtudes of Sij — > Sij + n"' and Sij — > + n"' derived from (3.4) are 



r- (3-8) 

with = — -^{ifj, + 'y^)75W. It is a simple matter to project out from (3.8) the strong 
decay amplitudes expressed in terms of component fields Bq, Bq, B^. For example, Eq.(3.8) 
leads to 

^ ^0 (3.9) 

A[{BQ)i^ ^ mijir^q)] = -%-u-MXl - X%)UQ {t < j), 

2v2/o 
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where use has been made of the identity (our convention being eoi23 = 1) 



■^i,v\K^^^^ = g^""!^ - g^^-f" + g^^-f^ - -f^Yl^, (3-10) 

and -^u^^ — u^, l^u^, — 0, v^u^ — 0. One-loop contributions are shown in Figs. 4 and 5 
respectively for S ^ Sn and 5" — > Tn. The resultant effective couphngs are given by 

{gi)es = gi ^ , (3.11) 

and 

[g2)eS = g2 7^, , (3.1ij 

^1 

where Z2{T) and Z2{S) are respectively the wave-function renormalization constants for T 
and 5" baryon fields, Zi and Z[ account for vertex renormahzation effects from Figs. 4a-4d 
and Figs. 5a-5d, respectively. 

To evaluate the renormalization constant Z2, we note that the self-energy amplitudes 
can be written as —i'E{k) and iTi{k)g^i, for T and 5" heavy baryon fields, respectively. Since 

E(^) ^ 8m- {Z^^ - l)v ■ k, (3.13) 



we find 



b -^0 

Z2mj) =i + lglJ2l [(^''^')- + + ^>^iA - 2A?,Ay 6,. (. ^ j), 

b 

for the antitriplet baryon with e^b^s being defined in (2.22). Likewise, we get 

,3 



(3.14) 



1 777 

^2(5y ) = 1 + 1 E i^akij + "^gkijW: 



(3.15) 



★ Our results for Z2{T) and Z2{S) [see also Eq.(3.18)] disagree with Eq.(3.3) of Ref.[17]. 
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for the sextet baryons Bq and Bq, where 



1 

t (3.16) 
Ctj =2 [(A'A^)n + (A^A^),, - 2A^,a5, - 2A^,a5,.(1 - <5,,)] . 

To work out the wave function renormahzation constant Z2 for heavy baryons, Eq.(2.32) 
and the following relations are useful: 

X^(AllAii)e^6 = ^(A22A22)M = ^tt + gCr/, 
b b 

XI ^'^13^31)^77' = ^iA3A2hn» = 2eK> (3.17) 

b b 

y^'^Ai2A2i)e7r'' = 2e7r, X(A33A33)e7r(> = -e^. 

It follows from Eqs.(3.14) and (3.15) that 

Z2(AQ)=l + 9^i(3e^ + e^), 

Z2(Sq) = 1 + ^^llSe^r + lOe, + 3e^), 
9 2 

Z2(Eq) = 1 + -^i(4e^ + 2e^ + -e^,) + 3^|(e^ + e^), (s^g) 

9 3 13 
M-q) = 1 + 85'i(2^'^ + + 6^'') + jgii^^T, + 2e^ + 3e^), 

9 o,. 8 



Z2(nQ) = 1 + -^f (4e^ + -e^) + egie^. 

As for the vertex diagram Fig. 4a for 5'^ — > Sij + n (since mass differences are generally 
small among different spin multiplets of heavy baryons, we will thus only consider pion 
emission due to the small available phase space), the general expression is 

Zi(4a) = 1 + E [lahj - ^9ki?) (3.19) 

with (see Appendix B for a derivation of the SU(3) group factors) 

Cij = g |(A^A^A^)jj + + X^ji{X^X^)ij + {X^X^)jiX'-j 

+ [Xl{X'X% + {X'^X')jjXi + a5,.(A^A^)u](1 - Sij) + (i - i)}, (3.20) 

Cij = ^{(AW),-,- - XUX'X% - Xl{X'X%{l - Sij) + (z - j)}- 
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Note that Z\ in Eq.(3.19) is worked out for 7r° emission, but due to SU(2) symmetry it should 
be also valid for charged pion emission. We next express the vertex corrections explicitly for 
5(SW) ^ 5(EQ)7f and 5(S'W) ^ ^(S'q)?: 

Zi(4a) = 1 + i^-g\ + 6^i) (e, + ^e,) + ^^^6, (3.21) 
for -5(E^*^) ^ -^(Eq)??, and 

Zi(4a) = 1 + ^ ^? (^-1% + 2e^ + ^e^^ + Qe^ + + ^e^^ (3.22) 
for 5(5'^^) ^ -S(5'Q)7f. Similarly for Sij Tij + n, Fig. 5a leads to 

Z[{5a) =1 + Y: {-^9U^j + ^9kl?) -^^^ {^ < j), (3-23) 



with 



Ci- = ^{(^'^'^')- - A^.-(A'A'')n + A^,(A3a''),-,- - (A^A3),,a5,. - (z ^ J 



Explicitly, 



(3.24) 



Z;(5a) = 1 + ^g! (^e, + ^e,^ + 3^i (^e, + ^e,^ (3.25) 



for ^(E^*^) ^r(AQ)7r and 

Z[{5a) = 1 + ^^1 Qe^r + + ^e,,) + 3^1 (^-^e^ + + ^6,^ (3.26) 

for 5(S'W) ^ T(SQ)7r. 

The cancellation between the wave-function renormalization of the Goldstone bosons 
and the 5-point vertex diagram for heavy-meson strong decays also persists in the baryon 
sector. This can be understood from the fact that the 5-point vertex in both cases arises from 
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expanding the chiral field Afj, to the third power in M and hence has the same structure. 
Therefore, Zi{Ad) — Z[(5d) — ^/Z2(n°-). Just hke their counterparts in the meson sector, 
all the seagull diagrams in Figs. 4 and 5 also do not contribute to the strong decays S ^ Sn 
and S Ttt to the lowest order in heavy quark expansion. First of all, we see from Eqs.(3.4) 
and (2.24) that the four-point vertices SSttH^ and TTtt" TT^ arc proportional to {qa — %) ■ v, 
where is the 4-momentum of the 71°". Second, there is no q dependence in the denominators 
of the amplitudes for Figs. 4c and 5c as the loop momenta can be chosen to be / and / + k 
for the light meson and heavy baryon, respectively. As for Figs. 4b and 5b, the q dependent 
terms in both numerator and denominator are always of the form {q ■ v). As a consequence, 
it will not give rise to the desired amplitudes shown in Eq.(3.8) after loop integration. We 
thus conclude that 

and 

{92U - 92 ^,^5^^ . (3.28) 

We now switch to the electromagnetic (Ml) decays 5" — > + 7 and -S" — > T + 7. The 
most general gauge invariant Lagrangian for magnetic transitions of heavy baryons given in 
paper II (II. 3. 9) can be recast in terms of superfields'*' 

£g) = - i3aitr(5^QF^^5'^) + Vla2t^^a(5ii{S^'Qv'' F'^^T) + h.c. 

3 _ (3-29) 
+ 3a'i tr(^^Q'(7 • FS^,) - -a'itr(rQV • FT), 

where a- F = a^^Fi^^ ^ and we have apphed Eqs.(II.3.47) and (II. 3. 61). The Lagrangian C-^' 
is also the most general chiral-invariant one provided that one makes the replacement 

Q^\{eQ^ + m^), Q'^Q'- (3.30) 

Note that contrary to Eq.(II.3.9), the Dirac magnetic moment terms do not appear in (3.29) 
because the Lagrangian is expressed in terms of velocity dependent fields and hence does 



★ The sign of the last two terms in Eq.(3.29) is opposite to those of the corresponding terms in Eq.(18) 
of Ref.[7]. 
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not contain Dirac magnetic moments to the lowest order. In the quark model the magnetic 
moments ai and 02 are simply related to the Dirac magnetic moments of the light quarks 
whereas a'^ is connected to those of heavy quarks. Explicitly, a'l is fixed by heavy quark 

symmetry to be ^^m^ ^^'^ ^^^^ only l/mg effect included for heavy baryon decays. In 

contrast, the couplings ai and 02 are independent of the heavy quark mass and spin. 

It follows from (3.29) that the tree amplitudes of heavy baryon radiative decays are given 

by 



A[S^.{v) ^ T,,- + 7(£, k)] = - 2^-a2 e^^apu-^v^k^e^iQii - Qjj)U^ (^ < j)- 

(3.31) 

For the radiative decay S ^ T + j, the two light quarks in the heavy hadron must undergo 
a spin-flip transition. Consequently, such decays will not receive any contributions from the 
magnetic moment of the heavy quark. As we will see later, this property persists at the 
one-loop level. The decay amplitudes in terms of component fields can be projected out 
from (3.31) by using the following relations: 



e^yai3Uv''k"'e^U^ ^iu{kn^ - e^^)j5U^ + -^ua^yk^e^'u, 

v3 

U^'ikyE^ - k^ei,)W ^vy{ki,e^ - k^ei,)u^ + -^u{k^^ - e^'lf)-i^u^ 

v3 

1 2 _ 

- -j=u'^{k^^ - £^^)75W + i-ua^yk'^e'^u, (3.32) 

2 

2 1 



To the lowest order, we have the following predictions from (3.31): 
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(3.33) 



2 



I 



where the superscript denotes the isospin quantum number. The above SU(3) relations will 
be violated in the presence of chiral-loop contributions as depicted in Figs. 6 and 7. Chiral 
corrections will modify the coupling constants oi, a'^ and 02 to 

(ai)eff =01—77^ + 5ai(6c+ 6d + 6e), 
^1 

7 ( Q\ 7 ('T\ 

(a2)eff - a2^^^^^|P^ + 5a2(7c + 7(i + 7e), 
^1 

where Zi (Z((5'5'7)), and Z[{TTj) are the vertex renormalization constants induced by 
Figs. 6a-6b, 7a- 7b and 8, respectively. In (3.34) we have anticipated the results [see (3.38) 
and (3.40) below] 

Z2iS) = Z[{SS-f), Z2{T) = Z[{TT-f). (3.35) 

We remark that since S'g^ — > Sg^7 is prohibited at tree level, its effective coupling (a2)efF 
is defined in a similar way as '^'q*^ '^q7- 

A[S{E'q^^*^) ^ r(SQ^)7] = -2^{a2)eSue^uapv''k''e^U'', (3.36) 

with 

(a2)eff = Sa2 (7a + 76 + 7c + 7d + 7e) . (3.37) 

Having calculated the renormalization constants Z2{S) and Z2{T) earlier, we now pro- 
ceed to evaluate Figs. 6a-6b, 7a- 7b and 8 for the vertex renormahzation constants Zi, Z[, Z'(. 
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Note that charge couphngs for the vertices SS^ and TT^ in Fig. 6a and SS^ in Fig. 7a are 
proportional to {v ■ e) and consequently cannot induce magnetic transitions. The contribu- 
tions to Zi due to the magnetic SS^ and ST^ couplings can be summarized as follows: 



z,(ssr, 6a) = 1 + ^ (^^gji.j - Q^^gmCij^ 
z[{ssr, 6a)^i + J2l (^^I'eX- + IglC'rj^ e,., 



Qii ~l~ Qjj 

(3.38) 



with 



Cij = I {(A^QA''),-,- - (A^A^),-,Q,, + 

= \ {(A'A''),i + (A'^A'^),, + 2A5,Af,. + 2X%xUl - 5,,)} , 
= {(A'A^).i + (A''A^)n - 2A5,A^^. - 2A5^.A^(1 - 5^)} , 



(3.39) 



for Sij Sij + 7, 

9 

b 

for Tij Tij +7,* and 



Z; {TTr, 8) = 1 + J ^ si [(A»A»)ii + (X'X^j + 2\lj\% - 2\i\''j,] ^ (i < j) (3.40) 



Z;'(.T,;7„)=I + i:l-2^,«.5^ + ^,i^j.,. ,<A (3.41) 



(3.42) 



with 

l^i = ^[ (A'QA"),, + (A''A''),,Q^^ + 2A^a5,Q,^ - ^ j) ], 

for Sij Tij + 7. It is easily seen from Eqs.(3.38), (3.40) and (3.14-3.15) that (3.35) is 
confirmed; the exact QCD result that does not get renormalized [13] is thus realized in 



★ Since the magnetic TT7 coupling vanishes in heavy quark limit and chiral corrections preserve heavy 
quark symmetry, we can be sure that no other diagrams such as the seagull diagrams and the diagram 
like Fig. 8 but with a photon attached to the light meson can contribute to Zi{T T7). We have 
checked this explicitly. 
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chiral perturbation theory. Note that when Qa — Qjj , the tree amphtude for 5" — > T + 7 
vanishes [see Eq.(3.31)] and Z" is undefined. In this case, one should apply Eqs.(3.36) and 
(3.37). The explicit expressions of the vertex renormalization constant Zi for each radiative 
decay mode are then given by 



ai 

4 4 5 ai 

ai 

, _i , 9^2/^. , , 9„_ 02^ 

ai 

3 3 ci 
^ ^. 6a) = 1 + -gfi-e, + 2e^) + 9^i^2 ^(e,), 

(3.43) 

and 



^i(Sq ^ 7; 6a) = 1 + g5i(^e7r + + g^,,) + 2^152 —{en + e,), 

^i(^g*^ 7; 6a) = 1 + 75i(e^ - ^e, + ^e^) + 9gig2 ^(e,), 

^i(^Q^^*^ ^ 7; 6a) = 1 + ^^i(^e^ + + ^e^) + 9^i^2 '^{^n), 

Q *^ ^ -Q 7; 6a) = 1 + i(-^e77 + + + ^5i^2^(-e7r + 2e^ + e^), 

Zi(5'q^(*) ^ S'q^ 7; 6a) = 1 + ^^??(7e, + ^e,) + ^^1^2 ^(e. + 6,), 



^i(^Q*^ ^ Ag7; 7a) = l + ^^i^2^(e7r + ^eK) + 3i/i(e7r + ^eK)> 

if*) „i „ . 9 ai,l 13 1,39, 

^1(4?' ^ ^q7; 7a) = 1 + -^i52-(2e7r + ye, + g^^) + i^K^Tr + 2e^ + 3e^), 

5a2(S'g^^*^ Eq^ 7; 7a) = ^^i^2ai(-e^ + e^) + ^gla2{-en + e^) 

(3.44) 

for 5" — > T + 7. Recall that the heavy quark magnetic moment a[ does not contribute to the 
5" T7 transitions at the tree level. Here, we see that this consequence of the heavy quark 
symmetry is preserved at the loop level. 

The electromagnetic five-point vertices in Figs. 6b and 7b are generated by the La- 
grangian (3.29) with the replacement Q ^ [M, [Q, M]] which comes from (3.30). The 
results are 

ZiiSij ^ Sij + 7; 66) = 1 + i 5] {[A^ [Q, X%i + (^ 

b 

b 
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i)} — - 



Oil ~l~ Qjj 



(3.45) 



■JJ 



Explicit expressions for Zi{6b) and Z'l{lh) for individual decay modes read 

°W^E?,7;66)= l-3e., (3.46) 











^ 7; 6&) 






- 7; 66) 





and 

^ 5| ^. 7^) ^ 1 _ + 26 J, (3.47) 

5a2(S'g'^*'* ^ 7; 76) = (-e^r + e^)a2. 

Figs. 6c and 7c arc the Feynman diagrams in which a photon couples to light pseudoscalar 
mesons. Care must be taken in order to ensure that the radiative decay amplitude is gauge 
invariant. As an example, we consider the decay Sg"^ ^Q^7 through the charged pion 
loop and the sextet baryon intermediate state. Its amplitude is given by 



where we have kept the superfield notation. We see that, aside from the k'^ term in the 
numerator, the {I ■ k) term arising from the first denominator in (3.48) also contributes; that 
is, they are all necessary for the sake of gauge invariance. By applying Eqs.(A3) and (A5) 
we arrive at 

MS; 6c) = z-^ !^ Wik^e^ - k^e,) W. (3.49) 
256 TT /q 

Likewise, the amplitude for the charged pion loop with antitriplet baryon intermediate state 
reads 

A,{T- 6c) = i A !^ ^ Tr{k,e^, - k^e,) W. (3.50) 

6Z TT /q 

Therefore, the radiative decay amplitude due to the sextet or antitriplet intermediate baryon 
state is separately gauge invariant, as it should be. Projecting out the spin-^ final state from 
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Eqs.(3.49) and (3.50) gives 



A,(EJ+1 ^ E+V; 6c) = ^ (^^? + - ^M^)75^^^- (3-51) 



Comparing this with the tree amphtude 



A(Eg^i ^ E+i7)tree = + 6a'ieQ)iZ(A;^9' - 6^^)751^^ (3.52) 



leads to 



<5ai(Ej+i ^ E+S; 6c) = ^(^,? + ,i)^ (3.53) 



for charged pion loops. On the contrary, the aforementioned {I ■ k) term does not contribute 
to the decay 5" — > T + 7 due to the presence of the totally antisymmetric tensor e^yajS- The 
general results including all charged meson loops are (see Appendix B for a derivation of the 
SU(3) group factors) 

= E al^^ (if?fe + grig-. 



h 



(3.54) 



'-33 



with 



-|<|([A', Q]A^^ + (^-J)^ 



Cy = -^{([A^ Q]X%i + {i^3)Y (3.55) 
el,- = - ^{([^'' + [A^ QlyAj, - (^ - j) 
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Working out the above general expressions, we obtain 

5ai(E+^W ^ S+S;6c) =5ai(S'-^^*^ S'^^S; 6c) 

o9l+ 92] I -?2" + ~F2" 



3e /3 2 2\ f^A ^^-^^^ 



and 



" 327r V8^^ ' V/o' 

5a,{^^ ^ ^ Aq 7; 7c) = - ^ ^-^ + - - 



'0 ^ ^0 



4 /n^ 



5a2(^Q ^^Q7,7c)--^(^--^ + -^j, (3.57) 

„-i „ s 3e gig2 ( 



^ --(*) r-,-i 

Note that 5a2 for the decay S'g^ — Sg^ is defined in Eq.(3.37) and it does not contain 
the factor of l/(Qi-j — Qjj)- We leave it to the reader to check that, in the SU(3) limit, all 
above results (3.43)-(3.57) do satisfy the SU(3) relations given by (3.33). 

It remains to see if the seagull diagrams are important. The relevant electromagnetic 
four-point vertices are described by 

--^ ^-^e^^p.S^^v^A^iQ, M\S- + ujl^-^S^A.iQ, M]T, (3.58) 
2V2 Jo V 2 /o 

where uses of Eqs.(3.4) and (2.49) have been made. Evidently, the vertices SSn"'^ and 
STn"'^ are independent of the photon's momentum k. This feature allows us to conclude 
immediately that Figs. 6e and 7e cannot cause magnetic transitions for the same reason as 
mentioned before for the case of 5" ^ S{T)'y. The remaining two seagull diagrams Figs. 6d 
and 7d also can be shown vanishing because of the identity w^m'^ = 0. Hence, 

Sai{6d + 6e) = Sa2{7d + 7e) = 0. (3.59) 
Finally, we note that the vertices SSn"' and STn"' in the Lagrangian (3.4) may suggest 
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that there is a possibihty of mixing between the sextet and antitriplet baryons through loop 
effects. However, this is not the case. An exphcit calculation shows that such loop diagrams 
vanish. The only possible mixing between S and T requires an e^jya^ tensor, but there are 
not enough variables to construct a Lorentz invariant. 

In summary, taking into account the mass differences among the Goldstone bosons, we 
found that SU(3) relations are generally broken in the strong and radiative decays of heavy 

hadrons. The leading chiral-loop corrections have nonanalytic dependence on niq of the form 

1 

rriq or niq Inm^. 
4. Applications 

In Sections 2 and 3 we have presented a general analysis of SU(3)-breaking effects in 
chiral perturbation theory for the strong and Ml radiative decays of heavy mesons and 
heavy baryons. In this section we will apply the results obtained so far to some selected 
decay modes of charmed mesons and baryons. Specifically, we choose the radiative decays 

D* D'j, Sc Ac7, — i> Sc7, S'* — S^7, S* — Sc7, ^c7 and the strong decays 

D* Dn, Sc — ^ AcTT as examples of applications. The decay rates of many of these modes 

have been explicitly calculated in papers I and II. Special attention is paid to the SU(3) 

1 

relations (2.45) and (3.33) for radiative decays. There are nonanalytic m| and mqlnniq 
chiral corrections which are responsible for SU(3) violation in chiral preturbation theory. In 
paper II, the relevant couplings are obtained from the quark model and they are related to 
the constituent quark masses. Hence, SU(3) violation is already incorporated there. In this 
section, a comparison between the two different approaches for SU(3) violation is made. 

We begin with the D* Dtt decay. When chiral loop corrections are included, we recall 
from Sec. 2 that its decay amplitude is given by* 

AiD* Djn-) = MyM^(A«),^.(£* . q), (4.1) 
JO 

where g^^ has been defined in Eq.(2.30) and q is the pion momentum. The decay widths 

★ Chiral logarithmic corrections to the decay D*~^ D^tt"*" are also calculated in Ref.[19] with the up and 
down quark masses being neglected. We agree with the results of Ref.[19] if the same approximation 
is made. 
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implied by the amplitude (4.1) are 



2 



(4.2) 



Note that D*^ -D^tt'^ is prohibited by isospin symmetry and D*^ — > D^n^ is kinemat- 
ically forbidden. The unrenormalized masses Mp and Mp. appearing in (2.1), (2.2), (4.1) 
and (4.2) can be inferred from the experimental measurement of heavy meson mass differ- 
ences. First, the SU(3)-invariant as term in Eq.(2.18) can be absorbed into a redefinition of 
Mp and Mp*. Then, (2.16) and (2.18) yield 

<7 -^o- M'^ll - M^. = -airrid - 7^ ( 3^ + 2:^^ + ^-p^ ) , (4.3a) 




and 



M^^ - = M^^^ - M,. = -arms - ^ U-^ + 3^^ ' " ^^"^^^ 



D7 " D' 



Eq.(4.3) imphes the mass difference^ 

phys phys phys phys (7^ I TTV^ TTl^ , 

M " - M [ = M \- M \ = -aUms - ma) - — i -3^ + 2^ + ^ • 4.4 



Note that formula (4.4) is different from that of Ref.[20] by a factor of 4/3 for the second term 

phys phys 

on the r.h.s. The parameter ai can be determined from the measured value M ^ — M^^ — 
99.5 ± 0.6 MeV [21], the current quark masses = 199 MeV, ma = 9.9 MeV [22], and the 
value of g chosen. The unrenormahzed masses and M^, can then be solved from (4.3).''" 



f In the presence of l/mq corrections, the heavy quark symmetry relations a2 = — ai and g = f/2 are 
no longer valid and Eq.(4.4) is modified to 

,,phy8 ,^phy8 ^ P ( n'^l ^"^K 



phys phys 1 /8 2 1 / ^'^l "'i) 



X At the preferred value g = 0.52 as shown below, we obtain = 1971 MeV and M^, = 2112 MeV. 
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We digress here to comment on the quantity 



^ (M M^J - (M -M^J= (M -M )- (M^,^ - M^^ ) (4.5) 

s s s s 

with the experimental value [21,23] 

= (0.9± 1.9)MeV, (4.6) 

which has recently received a lot of attention. Obviously, A^ vanishes in either the heavy 
quark or chiral limit. Beyond these limits, A^ receives two types of contributions. The first 
type is from the tree contribution induced by the l/rric corrections to the relation a\ — —a^- 
The second type comes from the self-energy one-loop diagrams with an insertion of lowest 
order mass splitting AM'^P{v)P'^ (v) and the ones without the above insertion but taking 
into account the splitting of g and //2. This was discussed in detail in our earlier paper [13] 
on heavy quark symmetry breaking effects. 

For the radiative decays D* — > D^, we see from Eq.(2.55) that the leading chiral correc- 
tions are dominated by Fig. 2e, and the wave function and vertex renormalizations. Before 
making concrete estimates on the parameter pi [cf. Eq.(2.41)], recall that in paper II we 
have considered the constituent quark model and incorporated SU(3) breaking into the hght 

quark charge matrix 

\ 



n -IMu J 

^ 3 / 



(4.7) 



To avoid any confusion with the current quark mass niq, we have used the capital letters to 
denote the constituent quark mass in (4.7). Therefore, the "effective" parameters {/3es)i are 
identical to and respectively for — > + 7 (i = 1,2,3) and hence in the 

quark model 

pr^ey<^(^ + f-l). (4.8) 
Using the constituent quark masses from the Particle Data Group [21] 

M^= 338MeV, Md= 322MeV, M^^SlOMeV, (4.9) 
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we get 



(4.10) 



for rric — 1.6 GeV. The effective coupling pi enters the radiative decay rate via 



with k being the photon momentum in the cm. frame. 

Recall that the quark model predictions agree very well with the two existing data on 
the D* decays [10]: the upper limit on the total width of D*'^ [12] and the branching 
ratio measurements by CLEO II [11]. It would be very important to see how well chiral 
perturbation theory does in this regard. There are four unknown parameters in the theory: 
g, P, rric and A, and available data do not permit an unambiguous determination of these 
parameters. Furthermore, other corrections such as l/mq effects may be of comparable 
importance as the chiral loop corrections. Unfortunately, these l/mg corrections are difficult 
to estimate [13]. Consequently, we will take a modest approach. We will set rric = 1.6 GeV 
and A = ~ 1 GeV, the chiral symmetry breaking scale. Then for (3 = 2.6 GeV^^ and 
3.0 GeV^\ wc let g take several different values at g = 0.5, 0.52, 0.60 and 0.75. The results 
are presented in Table I and Table 11. We see that for the total width of D*~^ to be of order 
130 keV or less we must have g ^ 0.52. On the other hand, the branching ratios for the 
D*^ measured by CLEO II appear to favor g — 0.52 and /3 — 2.6 GeV""*^. For this particular 
choice of g and /3, the predictions are quite close to the quark model results except for the 
radiative decay of D^'^. The decay rate for — > -0^7 is larger by an order of magnitude 
in chiral perturbation theory than in the quark model. Specifically, we find from Eq.(2.55) 
that 

pi = 4.59 e, p2=-1.14e, p3=-l-79e, (4.12) 

which should be compared with the quark model prediction (4.10). Of course, there is no 
reason to expect that these two different approaches for SU(3) violation should agree with 
each other exactly. Loosely speaking, SU(3) violation is treated nonperturbatively in the 



■k Strong decay widths are determined from (2.30) and (4.2) while radiative decay rates are calculated 
from (2.55) and (4.11). 
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nonrelativistic quark model, while it is calculated in terms of a perturbative expansion in 
chiral effective Lagrangian theory. Nevertheless, the fact that the predicted decay rates and 
branching ratios for D*^ and agree so well with the quark model and the data suggests 
a consistency between "model and theory" . It will be extremely interesting to measure the 
rate for L>*+ D^l to see which prediction, if any, is closer to the truth. It is obvious 
from Table I that a smaller total width of the D*^ [12] can be easily accommodated without 
upsetting the agreement of the predicted branching ratios for both D*^ and D*^ with the 
CLEO II data [11]. 

At this point it is worthwhile to compare our favored choice oi g — 0.52 and (5 — 
2.6 GeV""*^ in chiral perturbation theory with the corresponding parameters g^^y^ and I3^y^ in 
the quark model. Recall that g^^^^ = 0.75 [3] and (3^-^^ = l/M.^ = 2.96 GeV""*^ [10]. Since g^-^j^^ 
is nonperturbative in nature while g is an unrcnormalized parameter in the chiral Lagrangian 
approach, they are a priori not the same. To see their relation, we write down the D* — > Dtt 
amplitude in the quark model: 



A(D* ^ DjTT'') = h^JM^'^'^'M'^r Wji^* ■ q)- (4.13) 

For this to be in accordance with (4.1) predicted by chiral perturbation theory, we are led 
to 



Since Z2 > Zi, fn > fo, > ^ and M^, > M^/ , it is evident that g < g^-^^. By the 
same token, the relation between /3 and /3q^ can be understood along the same line. 

We next shift our attention to the heavy baryon sector. For the strong decay of heavy 

baryons, we shall only consider the decay Ec — > AqTT which is experimentally seen, although 

f*) 

its rate has not been measured. (It is likely that no other ' B^tt or Bq — > Bqtt strong 
decays are kinematically allowed.) The decay rate for — > AcTT is given by [3]^ 

^A^TT)-— q, (4.15) 



t Note that fo is replaced by /,r in Ref.[3], since there we work in the tree approximation only. 
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where [see Eq.(3.28)] 



which includes chiral loop effects. Before examining the SU(3)-breaking effects induced by 
chiral loops in the electromagnetic decays of heavy baryons, it is instructive to examine the 
tree amplitudes of (3.33) to see what SU(3) violations arc expected from the constituent 
quark model. From Eqs.(3.31) and (4.7) we find the following quark-model predictions: 

(a.)S'(sC ^ Etj) =aT + OMaT = O.SSeGeV"', (4.17) 

(«2)eff c ^ -cV = 3«2 ( Ms ~ ) = -0-13«2 = -0.078 eGeV , 

and 

(ai - 8a'i)eff (E*° ^ EO7) = — ^ - 8ai = 1.05ai - 8a[ = -1.45 eGeV~\ 



/ ^ QM , ^ , QM M^i , QM , 1 

(ai - 8ai)cfj (f]* ^ Qcl) =^1 - Sa'i = 0.66ai - 8a[ = -1.07eGeV"\ 
(ai - 8ai)eff (^ , - ^ ,7) = + mT j " 

QM , 1 

= 0.86ai - 8ai = -1.26eGeV"\ 
(«i + 16oi)eff (= c ^ = c7)=3ai (g-g^^J+l^ai 

= 1.34ai^' + 16a'i = -0.49eGeV"\ 

QM , 1 

= 0.95ai + 16a'i = -0.104 eGeV"\ 

where we have made use of the quark model predictions [10] 

QM e QM e 

with = 1/-/W^- The parameter a'j^ is predicted by heavy quark symmetry to be [10] 
Prom the lowest order Lagrangian we have the SU(3) relations of (3.33). It is thus clear that 



(4.18) 
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a measurement of the decay rates of 

[V ^ ^c7, ^-c7, i^c^^^cT, 

(ii) H'° ^ H°7 relative to S'^ ^ S+7 or E+ ^ A+7, 

(iii) S':+ ^ S'+7, S*+ ^ E+7 
will provide a nice test on SU(3) breaking. 

We are now ready to compare the quark model predictions with those from chiral pertur- 
bation theory using (3.27-3.28) and (3.34) for heavy baryon's decays. There are six unknown 
parameters in the theory: A, nic, gi, g2, ai and 02. Unfortunately, there does not exist any 

experimental data for the heavy baryons to guide the choice of these parameters. It is clearly 
premature to make an all-out effort to study how the predictions vary with the parameters. 
To gain a general idea of how chiral perturbation theory is compared with the quark model, 
we minimize the number of parameters by setting A = 1 GeV and rric = 1.6 GeV as before, 
and by making use of the quark model relations [10] 

1 [2 e e 

51 = 92^ -y ^9, ai = a2=^/3. (4.21) 

The remaining two parameters g and (3 are finally fixed with the values favored by the 
available D* decay data: 

£, = 0.52, /? = 2.6GeV"^ (4.22) 

The relevant wave function and vertex renormalization constants Z2, Z[ and Z" for various 
processes are given by (3.18), (3.25-3.26) and (3.41-3.47). Numerically, 



(4.23a) 



Z2(Ac)=1.40, Z2(5c)=1.83, 

Z2(Ee) = 1.12 , Z2(S',) = 1.14 , Z2(Qc) = 1-19 , 
and 

Zi(Ec ^ AcTt; 5a) = 1.07 . (4.236) 
From (4.15), (4.16) and (4.23) we find that 

r(E++ ^ A+7r+) = 1.89 MeV (4.24) 

in the chiral Lagrangian approach as compared with 2.45 MeV in the quark model [3]. 
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Vertex corrections to radiative decays are 



Zi(S*++ ^ S++7) =0.77, ^ S+7) = 0.62, ^ £^7) = 0.92 

Zi(S':+ ^ 5'+7) = 1.00 , Zi(H':° ^ H'°7) = 0.78 , ^ Oe7) = 0.63 , (4.25) 

Z'iiK^ - Ae+7) = 0.88 , z;'(S'r ^ S+7) = 0.87 , 



and 



and 



(ai + 16ai)eff(H'*^ ^ 5';^7) = 1.02ai + 16a'i = -0.05eGeV-\ 

(ai + 16ai)eff(S*+ ^ S+7) = 1.35ai + IQa-^ = -0.35eGeV~^ 



(4.26) 



5ai(Ej^^*^ ^ S+S;6c) = 5ai(H'g^^*^ ^ S'^S; 6c) = 0.24eGeV-\ 
5ai(r][j^ ^^Q7;6c) = 5ai(SQ^*^ ^ S^7;6c) = 0.38eGeV-\ 
(5ai(Eg^^*^ ^ Sq^7;6c) = (5ai(H'|^*^ ^ 5'|7;6c) = 0.11eGeV"\ 
5a2(Eg^*^ ^ Aq7;7c) = -0.08eGeV"\ 
5a2(S'|^*^ ^ s|7;7c) = -0.16eGeV~\ 
^02(2'^^^*^ ^ Sq^ 7; 7a + 76 + 7c) = 0.15a2 - 0.03 eGeV"^ 

We get from Eqs.(3.34), (4.23a) and (4.25-4.26) that 

(ai - 8a'i)eff(S*° ^ E°7) = 1.08ai - 8a'i = -1.36eGeV"\ 

(ai - 8ai)efr(l^c ^cl) = IMai - Sa'i = -1.68 eGeV"\ (4.27a) 
(ai - 8a'i)eff(5'*° ^ 5'c7) = - 8ai = -1.45eGeV"\ 



(4.276) 



A comparison between (4.27) and (4.18) indicates that the two approaches give similar 
results for some processes such as E*o ^ EO7 and S'*° ^ S'°7, while they give dramatically 
different results for others such as E*+ — > E+7 and S'*"*" — > S'+7. For the decay amphtude 

^1(^6* ^ + 7) = my'ik^,^ - £M^)75^i^, (4.28) 
with r)2 being given by Eq.(3.33), the corresponding decay width is [10] 

k ( \ ^ 

V{Bl - ^6 + 7) - ^ (1 - ^1 (3mf + mj), (4.29) 
where mj (m/) is the mass of the initial (final) baryon in the decay. This formula is presented 
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here for completeness, and we will not give any predictions for the decay rates of spin-| heavy 
baryons as their masses are still unknown. 

As for the electromagnetic decays Bq ^ + 7, Eqs.(3.34), (3.44) and (4.23-4.26) lead 

to 

(a2)efF(S+ ^ A+7) =1.2702 = 0.67eGeV-\ 

{a2)es{^'t ^h) = l-36a2 = 0.72 eGeV-\ (4.30) 
(a2)eff(5'c ^ S°7) =0.10a2 = 0.05 eGeV-^ 

Note that 5'^ 5^7 receives two contributions: one from Fig. 7a and the other from Fig. 7c 
and there is a cancellation between the two. The sign of this amphtude in (4.30) is actually 
opposite to that in the quark model [cf. Eq.(4.17)]. With the formula 



r(56^i^3 + 7)= -v^k^ (4.31) 

TT 



for the radiative decay amplitude 

A{B6 ^ S3 + 7) = r]lu-sa^„k^'e''u6, (4.32) 
the decay rates of various decays are then predicted to be* 



r(s+ - 


- Ac+7) 


= 112keV, 


r(s': - 


- S^7) 


= 29keV, 


r(s'° - 


-SO7) 


= 0.15keV, 



(4.33) 



which should be compared with the corresponding quark-model results: 93 keV, 16 keV and 
0.3 keV. 

The predicted SU(3) breaking patterns in chiral perturbation theory for radiative decays 



★ All the charmed baryon masses are taken from the Particle Data Group [21]. For the mass of ^'c, we 
employ the hyperfine mass splitting — ~ 100 MeV derived in Ref.[25]. 
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of charmed mesons and baryons are [see (4.12) and (4.27)]: 



pol pol 



pol pol pol 



2 



(4.34) 



|2 



E i^(sr - E+7)r » E |^(-'c^ - -'c^), 

pol pol 

which are opposite to what expected from the quark model [see (4.10) and (4.17-4.18)]. It is 
thus extremely important to measure the decay rates of the above-mentioned modes to test 
the underlying mechanism of SU(3) violation. 

We should reiterate the tentative nature of the discussion in this section. In the heavy 
meson case, the data for D* are useful but not decisive in selecting the values of the param- 
eters. Other important effects such as 1/mq corrections have not been incorporated. The 
procedure employed here is by no means a best fit. The apparent agreement between the 
quark model and chiral perturbation theory is pleasing, but it should be taken with caution. 
In the heavy baryon sector, the situation is even more difficult. First of all, we do not have 
any experimental data. Secondly, there are many more parameters here than in the heavy 
meson sector. The choice of parameters based on (4.19)-(4.21) is necessarily ad hoc. The 
only virtue is that it provides a rule to fix the parameters. It will not be surprising if nature 
picks a very different choice from the present one. Thirdly, the renormalization effects are 
very significant, as evident from (4.23) and (4.24). We may question the validity of chiral 
perturbation expansion.^ Fourthly, large cancellations occur sometimes. For all these rea- 
sons, the results obtained in this section should not be considered final. Rather, they make 
it clear that more experimental and theoretical works are needed. 

5. Conclusions 

As we have emphasized before, the lowest-order low-energy dynamics of heavy mesons 
and heavy baryons is completely determined by the heavy quark symmetry and chiral sym- 
metry, supplemented by the quark model. In order to gain a different perspective of our 



t The similar problem in the meson sector is less severe except for the D*^ decay. For g = 0.52, we have 
Z2{D*°) = Z2{D*+) = 1.18, Z2{D*+) = 1.32 and all the Zi's are very close to unity for the strong 
decays, and Zi(£)*° ^ £)°7) = 0.74, Zi{D*+ D+^) = 0.87, Zi{D*+ D+-f) = 0.57 for radiative 
decays. 
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previous quark model calculations, we have investigated in the present paper chiral symme- 
try violation induced by the light current quark masses. Chiral corrections to the strong and 
electromagnetic decays of heavy hadrons are computed at the one loop order. The leading 
chiral- loop effects are nonanalytic in the forms of m| and ruq Irnriq. We have also confirmed 
in chiral perturbation theory at the one loop level the exact result in QCD that the heavy 
quark contributions to the Ml transitions of heavy hadrons are not modified by light quark 
dynamics [see (2.51) and (3.35)]. Moreover, all consequences of heavy quark symmetry are 
not affected by chiral corrections. We then apphed our results to the radiative and strong 
decays of charmed mesons and charmed baryons. 

We found in the meson sector that for a particular set of the parameters g and /? inferred 
from the measured branching ratios of the D*^ and D*"*" and the upper limit on the D*~^ 
rate, the predictions of the strong and electromagnetic decays of the D* mesons in the chiral 
Lagrangian approach are not very different from the quark model results except for the 
radiative decay of the Dg"*". Using the quark model to relate the four unknown parameters 
in the baryon sector to g and P, we have computed the strong decay Sc AcTt and the 
radiative decays Bq ^ Bq + and Bq ^ B^ 7. We found that the chiral-Lagrangian and 
quark-model approaches in general give similar results for many processes, whereas they yield 
drastically different predictions for others such as E*+ and 5'*^ Moreover, 

the predicted SU(3)-breaking patterns as shown in (4.34) are opposite to the quark model 
expectations. It is thus of great importance to measure the decay modes listed in (4.34) to 
test the underlying mechanism of SU(3) violation. 

A more meaningful and improved comparison between theory and experiment has to 
wait until more data on heavy hadrons become available. In the meantime, it is important 
to have a better theoretical understanding of the corrections due to the l/mg effects. We 
believe that both SU(3) symmetry breaking and the l/niQ effects should be considered 
simultaneously. We will return to this effort sometime in the future. 
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APPENDIX A 

In this Appendix we list some of the useful Feynman integrals which are relevant to this 
paper: 



/ 



(27r)^ (^2 -m^ + ie){vl + ie) '"''^ IGn^ ^ m^' ^ ' 



I 



/ 



"^^1^ ^^m2{^^'"'' ~ '^v^v^){v ■ k + 5m)m^ - 2{g^,, - AvnV„){v ■ ~k)5w? I ^^^^ 

+ ^T^^9nv - 'ivnVv){v ■ k)5m, 
Idtt 



(2^ [(/ + A:)2 - m2 + ^e] - + ^e){y ■ I + re) (with = and ^ = 0) 

Tfi I 11 

= «^ S [(s'm!^ - VnV^)ka + cyclic in i/, a] + ^^k^k^ka > {A5) 
^327r2 ^ ^ 1 ^^^'"'^ ^^^(J'^'^J^a + cyclic m i/, a J - 2v^v„Va > , 



where the common factor 



1111, , . . 



with e = 4 — n has been lumped into the logarithmic term In A^. as A is an arbitrary 
renormalization scale occurred in the dimensional regularization approach. Note that only 
the leading contributions hnear in are retained in (A4). Our results for the integrals (A3) 
and (A4) are in agreement with Ref.[17]. Eq. (A5) is derived by assuming v-k — 0, where k is 
the momentum of the outgoing photon which couples to heavy hadrons. This is a legitimate 
assumption as long as both incoming and outgoing hadrons are on shell. To evaluate above 
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integrals, we first apply tensor decomposition as well as momentum expansions to reduce 
them into the following prototypes: 

. , [A^ 1 . . 

^^""^ J (27r)^(/2-m2 + ^e)«(^;-/ + ^e)' ^ ^ 

and 

T( ] = I (AO^ 

(27r)"(Z2-m2 + ie)«(T;-Z + ie)2- 

We observe that all the external momenta have disappeared from the denominator due to 
the expansion: 

1 1 /"^ 2l-k + k^\ 



(/ + A;)2 — 777.2 ^2 _ ^2 \^ ^2 _ .^2 ^ 

This expansion is a valid procedure since the momentum k is very soft in the current context. 
Similar type of expansions can be applied to Eq. (A4) for bringing up the residual momentum 
k to the numerator. 

To carry out integrals Ii and I2, the usual Feynman parametrization 

J [ax + 6(1 - ^ ' 



1 

1 _ r(m + n) r x'^-'^ii-xy-'^dx 

~ r(m)r(n) J [ax + 6(1 - 




does not allow one to shift the loop momentum I without destroying the structure of the 
integral. Hence, it is convenient to combine denominators using the identity 



00 



1 _ r(m + n) f 2"A"-^dA 

r(m)r(n) J (a + 2A6)"*+"' 



(^11) 



which is obtained from (AlO) by changing variable x — j^2X- Sometimes the integration can 
even be done without combining denominators. We consider the integral Ji as an example 
to illustrate the importance of the ie terms presented in the denominator. We write 

ill a 



(27r)" (/2 - m2 + ie)«(v ■ I + ie) 

1 fA^d-l 1/1 ^ ^ ^^^^^ 

+ 



2 J (27r)'^(Z2-m2 + ie)« V'y-^ + ^e -v ■ I + ie 
The second line follows from the fact that h is invariant under the substitution I ^ —I. 
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Since 



-iirSivl), {A13) 



V ■ I + ie V ■ I 
the integral becomes 

Since / is a scalar integral, we can easily evaluate it by going to the rest frame of v. In the 
rest frame of f , = (1, 0), we find 

h{a) = -m(-l)-A^ / '^-^ (A15) 

where I is a vector in the (n — l)-dimensional space. Performing the angular integration gives 



- m-2 



I^(a) = \ — { J- / dl-^ M16 



Now the /-integration is standard, which gives 



It is easy to see that I\{a) is finite for any integer a. Taking n = 4, we have 



To evaluate /2(Qi), it is more convenient to apply the identity (All) since the principal part 
of the integral no longer cancels. By the power counting argument, one expects /2(1) to be 
logarithmically divergent. However, as we will show momentarily, such ultraviolet divergence 
will appear in the A integration rather than in the momentum one. Therefore some care is 
needed in order to consistently treat the infinities. To illustrate this, let us first apply Eq. 
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(All) to combine denominators of hi^)- This procedure gives: 

oo 



72(1) = J SXdX J A' 



(27r)«[/2- (m2 + A2)]3- 



(A19) 



It is obvious that no infinities arise in the momentum integration. Carrying out the momen- 
tum integration, we obtain 



(^20) 



where we have changed the variable such that k — X"^. Apparently the ultraviolet divergence 
now appears in the k integration. To consistently implement the dimensional regularization, 
let us perform an analytic continuation on /2(1) with the aid of 

/ "^"^4^^^'^^^^^^ ^ / "^"^ ■ ^^wj^^^ ^""J ^^^^^ 

where / is an arbitrary function of I. The l.h.s. of Eq. (A21) can be set to zero since 
the integrand is a total derivative. Now -^2(1) can be defined in terms of ^2(2) through Eq. 
(A21): 

2m2 



^2(1) 



{A22) 



A-n 

Since /2(2) is convergent, the infinity in /2(1) manifests itself as a simple pole at n = 4. To 
determine /2(1) including its finite parts, one has to evaluate /2(2) to C(4 — n). Applying 
Eq. (All) gives 

00 



h{2) 



24:XdX 



1 







(27r)"[Z2-(m2 + A2)]4- 



{A23) 



Carrying out both integrations, we arrive at 

2i 



H2) 



By Eq. (A21), we obtain 



167r2m2 



^2(1) = 



4i 



167r2 



l-|(7.-ln.-21n2) + ilnA_ 



;(7,-ln7r-21n2) + lln^ 



(A24) 



(A25) 



To compute any hia) with a > 2, one may repeatedly apply Eq. (A14). The general result 
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is given by 



with a> 2. 



Ha) = \ (A26) 



Finally, the integral (A5) can also be carried out by combining the first two denominators 
via the conventional Feynman parametrization (AlO), followed by another combination with 
the third denominator {y ■ I + ie) using the identity (All). As a consequence, 



oo 1 

A I d\ I dx 



[{I + kY — rn^ + — rn? + ie){v ■ I + ie) 



1 



{A21) 



X 



[{I + + kx)"^ + x{l — x)k'^ — — m^]^ 
After the momentum integration, the integration over x becomes trivial. 



APPENDIX B 

In this Appendix we explain the procedures used in Section 3 to obtain the results 
for strong and electromagnetic vertices, especially how the SU(3) group factors are defined 
and arrived at. We will consider two examples: one for a strong vertex and one for an 
electromagnetic vertex. For the strong vertex, we will take Fig. 4a with sextet baryons as 
the intermediate states. The corresponding tree amplitude can be written as 

^i^ij ^ ^ki + = 4 J^'mfir^U''v''<l^W{SuMBh'^B) l^y), {Bl) 

where B is the sextet baryon matrix as defined in Eq.(3.2); \Sij) and 15";^^) are respectively 
the SU(3) flavor wave functions for the initial and final baryons, such that 

Bmn\Sij) — ~~7={5mi^ni + 5mj^ni)\^) {i ¥^ j)i / t> \ 

Bmn I ^ii ) — ^mi ^ni 1 0) ■ 

Hence for the specific case of Sij —>■ Sij + tt^, we have 

The loop contribution from Fig. 4a with sextet intermediate states can be written in the 
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following form 



b 

where M^y{h) is the contribution from the loop excluding the SU(3) matrices: 



m 



where we have neglected the residual momentum k of the initial baryon and the external pion 
momentum q in the propagators since we are only interested in the leading order contribution. 
Using the result (A3) for the integral, we obtain 

M^,(6) = -2ie^af3.v''q''j- , (^6) 

where e^b is defined in (2.22). The operator G{b) is a product of traces involving SU(3) 
matrices and the sextet baryon matrix B: 

G{b) = tT{Bh^B) ■ triBh^B) ■ tv{B^X^B). {B7) 

Note that the adjacent intermediate S's and in (B7) combine in pairs to form sextet 
"propagators" in the SU(3) flavor space, namely 

Thus G{b) can be rewritten as 

G{b) = ^tr ^B^X^'X'^BX"^ + B^X^X'^BX^^ + B^X''B{X''XY + B^X^X^'X^B^ , {B9) 

where the superscript T signifies the transposition of a matrix. Again, for the case of 
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Sij — > Sij + TT^, the flavor matrix element can be easily worked out to be 



{Sij\G{b)\Sij) 



{BIO) 



where ^ij is given by (3.20). When (B3) and (BIO) are combined, we obtain 

AsiS^j ^ S^.+n') = -2ij-e,af3.U''v%^U'' (^^i)' E^-^^^ " (^H) 
Adding up the tree amplitude (B3), we find 



A{Si:j^S^.+7:^) 



■ egi 



1 y 2Y- 

1 - 1^1 



■ (512) 



The quantity in the square brackets of (B12) is the contribution to Z^^ ^ due to the interme- 
diate sextet baryons of Fig. 4a. It agrees with the term proportional to g\ of Eq.(3.19). 

The contribution to Zi due to antitriplet baryon intermediate states of Fig. 4a can be 
computed similarly. The SU(3) fiavor "propagator" for an antitriplet is 

The difference in normalization between (B13) and (B8) is a result of how the matrices for a 
sextet and an antitriplet are defined [see (3.1) and (3.2) ]. The minus sign in (B13) accounts 
for the antisymmetry of the antitriplet baryon matrix T (or B^). 

We now give an example of the calculation for a loop diagram in which a photon couples 

to Goldstone bosons. Consider the Feynman diagram of Fig. 6c with intermediate sextet 
baryons. The interactions of the Goldstone bosons with the electromagnetic vector potential 
yl^ are introduced into Eq.(2.20) by the substitution 



a^E ^ D^E = a^E + ieA^[Q , E] 



(514) 



The electromagnetic vertex for tt** — > tt'' + 7 is given by 



Ain"" ^ TT^ + 7) = -ie{2q - k) ■ e 



--tr(A«[Q , A^]) 



{Bib) 



where q and k are the momentum of tt** and the photon, respectively. The quantity in the 
square brackets is equal to one when the Goldstone boson has a unit positive charge. It will 
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be part of the group factor defined below. Fig. 6c with intermediate sextet baryons gives a 
contribution to S ^ S + ^: 

MSrj - S^, + 7) = YJ^'M'^M U^{Su\G'{h)\Si^), (S16) 

h 

in a notation similar to (B4). We find 



2 ^ 



Notice that the Goldstone bosons tt" and tt^ on both sides of the electromagnetic vertex must 
have the same mass which we denote as uiT^b. The integrals needed in (B17) are given in 
Appendix A. After neglecting a contribution to the convection current proportional io v ■ e, 
we get the gauge invariant result 

= '^iik.e^-k.e,) . (518) 

The group factor operator G'{h) is given by 

G'{h) ^Yir{B^\^B) ■ tr(StA'*S) -^tr (^A«[Q , A^]) . (S19) 

a 

Making use of the sextet propagator (B8) and the identity 

^'jkKnn = '^{^jnSkm " -^^jk^mn) , (-B20) 

a 

we obtain 

G'{b) = -^tr |stA^[Q , X^]B + B^\^B[Q , xY] ■ (521) 
For the case Sij — >• Sij + 7, we find 

{S^j\G'{b)\S^j) =-l{ (x'^iQ , X']) + xUQ , X%{1 - 6ij) 
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As a consequence of the fact that Q is diagonal, we have 



{B23a) 



{B23b) 



Finally, 



{Sij\G'{b)\Sij) = ^ij , 



(S24) 



where is defined by (3.55). Combining (B18) and (B24), we obtain 



X 



m^^b 3 2f 1 



-72- ■ 7fi^y- 



{B25) 



The quantity in the square bracket gives rise to the term proportional to gf in (3.54) 
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Figure Captions 



Fig. 1. Chiral-loop diagrams contributing to the strong decay P* Ptt". For simplicity, the 

wave-function rcnormaUzation and mass countcrtcrms are not shown in Figs. 1-8, but 
the necessary mass renormahzation is to be understood. 

Fig. 2. Chiral-loop diagrams contributing to the radiative decay P* — > P7. 

Fig. 3. Amplification of Fig. 2e with all charged meson loops specified. 

Fig. 4. Chiral-loop diagrams contributing to the strong decay S Sn"'. In 4(a), the external 
pion may originate from an SStt or STir vertex. 

Fig. 5. Chiral-loop diagrams contributing to the strong decay S Tir"". In 5(a), the external 
pion may originate from an SSir or STir vertex. 

Fig. 6. Chiral-loop diagrams contributing to the radiative decay -S" — > S^. In 6(a), the photon 
may originate from an SS^ or ST^, or TT^ vertex with a coupling constant ai or 02 
or a'l- In 6(c), (d) and (e), the intermediate state can be an 5" or T baryon. A diagram 
similar to 6(b) but with a photon attached to the light meson does not contribute to 
the Ml transition oi S ^ S^. 

Fig. 7. Chiral-loop diagrams contributing to the radiative decay S Tj. In 7(a), the photon 
may originate from an SS^ or ST^ vertex with a coupling constant ai or 02. In 7(c), 
(d) and (e), the intermediate state is an S baryon. 

Fig. 8. Chiral-loop diagram contributing to the radiative decay T — > T7. This is the only 
loop diagram contributing to the Ml transition T T'-f. All other possible diagrams 
do not contribute; see the footnote after Eq.(3.40) for discussion. 
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Table Captions 



Table I. Predicted decay rates (in units of keV) and branching ratios (in parentheses) of strong 
and electromagnetic decays of charmed mesons for various values of g with j3 = 
2.6 GeV^^ and rric = 1.6 GeV. For comparison, the quark model predictions [10] and 
the experimental branching ratios measured by CLEO 11 [11] are given in the last 
two columns. The numbers under "quark model" differ somewhat from those given in 
Ref.[10] because of the more precise pion masses used here. 

Table II. Same as Table I except for /3 = 3.0 GeV"^ 
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Table I. Predicted decay rates (in units of keV) and branching ratios (in parentheses) of strong and el 
of charmed mesons for various values of g with /3 — 2.6 GeV""*^ and rric — 1.6 GeV. For comparison, the qm 
[10] and the experimental branching ratios measured by CLEO II [11] are given in the last two columns, 
"quark model" differ somewhat from those given in Ref. [10] because of the more precise pion masses used 



Reaction 



g = 0.52 g = 0.6 



g = 0.75 



quark model CLE 



D*+ 
D*+ 
D*+ 

D*+ 



D 



*o 



total 



total 



D*+ D+ 



7 



78.8 (68%) 87.9 (68%) 133.4 (68%) 276.0 (68%) 



35.7 (31%) 39.8 (31%) 60.5 (31%) 

1.9 (1.7%) 2.0 (1.5%) 2.1 (1.1%) 

116.4 129.7 196.0 

54.1 (61%) 60.3 (65%) 91.5 (76%) 

34.0 (39%) 33.0 (35%) 28.9 (24%) 

88.1 93.3 120.4 



4.6 



4.5 



3.7 



125.0 (31%) 
2.5 (0.6%) 
403.5 

189.3 (90%) 
19.9 (10%) 
209.2 

2.1 



102 (68%) 
46 (31%) 
2 (1.3%) 
150 

70 (67%) 
34 (33%) 
104 

0.3 
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Table II. Same as Table I except for (3 = 3.0 GeV"^ 



Reaction 



^ = 0.5 



g = 0.52 g = 0.6 



g = 0.75 



quark model CLE 



D*+ 
D*+ 
D*+ 
D*+ 



total 



D*^ total 
D*+ ^ D+i 



78.8 (67%) 
35.7 (30%) 
3.4 (3%) 
117.9 



8.3 



87.9 (67%) 133.4 (67%) 
39.8 (30%) 60.5 (31%) 



3.4 (3%) 
131.1 



8.2 



3.8 (2%) 
197.7 



54.1 (53%) 60.3 (56%) 91.5 (67%) 
47.5 (47%) 46.6 (44%) 42.6 (33%) 
101.6 106.9 



134.1 



7.5 



276.0 (68%) 
125.0 (31%) 
4.7 (1%) 
405.7 

189.3 (85%) 
33.3 (15%) 
222.6 

6.0 



102 (68%) 
46 (31%) 
2 (1.3%) 
150 

70 (67%) 
34 (33%) 
104 

0.3 
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